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SUMMARY 


This  report  summarizes  the  first  year  work  under  the  general  heading  of 
"Effect  of  Local  Material  Imperfections  on  the  Buckling  Behavior  of 
Composite  Structural  Elements". 

It  describes,  through  the  introduction ,  two  important  areas:  one  of 
delamination  buckling  and  growth  and  one  of  the  holes  and  foreign 
inclusions. 

In  the  subsequent  sections,  the  important  subject  of  delamination 
buckling  and  delamination  growth  is  presented  with  sufficient  detail.  The 
geometries  considered  are  flat  laminates  and  the  emphasis  is  to  establish 
the  load  carrying  capability  (damage  tolerance)  of  the  laminate.  It  is 
clearly  shown  that  two  important  parameter  govern  the  behavior  of  the 
delaminated  (damaged)  laminate:  the  size  of  the  delamination  and  its 
position  (especially  relative  to  its  distance  from  the  surface).  Depending 
on  these  two  parameters  the  damage  tolerance  of  the  laminate  is  either 
governed  by  buckling  or  by  the  fracture  toughness  of  the  material 
(delamination  growth).  Finally,  the  second  year  effort  is  described.  This 
includes  the  analysis  of  delaminated  cylindrical  shells  under  various  load 
conditions. 
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I.  INTRODUCTION 

The  constant  demand  for  lighter  and  more  efficient  structural 
configurations  has  led  the  structural  engineer  to  the  use  of  new  man-made 
materials.  At  the  same  time  this  demand  forced  upon  him  very  sophisticated 
methods  of  testing,  analysis  and  design,  as  well  as  of  fabrication  and 
manufacturing.  The  recent  explosive  progress  in  producing  and  using 
composite  materials  has  pointed  towards  the  clear  possibility  of  man 
creating  specific  materials  for  specific  missions.  At  the  same  time  it  has 
been  realized  that  some  severe  needs  are  created  (related  to  this  effort) 
namely  (i)  the  complete  understanding  of  the  behavior  of  composite 
materials  and  what  influences  this  behavior,  (ii)  the  establishment  of 
design  criteria  upon  which  proper  use  of  composites  rests,  and  (iii)  the 
training  of  engineers  in  the  area  of  composite  structures.  In  all  three 
items,  it  is  important  to  recognize  that,  with  the  advent  of  composite 
media,  certain  new  material  imperfections,  can  be  found  in  composite 
structures  as  well  as  the  better  known  imperfections,  that  one  finds  in 
metallic  structures.  Thus,  broken  fibers,  delaminated  regions,  cracks  in 
the  matrix  material  as  well  as  holes,  foreign  inclusions  and  small  voids 
constitute  material  and  structural  imperfections  that  do  exist  in  composite 
structures. 

I . 1  Delamination  Buckling  and  Growth  (Flat  Plates) 

One  of  the  most  common  failure  modes  in  composite  materials  is 
delamination. 

Delamination  is  developed  as  a  result  of  imperfections  in  the 
production  technology  or  due  to  the  effect  of  certain  factors  during  the 
operational  life  of  the  laminate,  such  as  impact  by  foreign  objects.  The 
presence  of  delamination  in  a  laminated  composite  material  may  cause  local 


1 


4 


« 


1C 


buckling  and  reduction  in  the  over-all  stiffness  of  the  structure,  which 
may  lead  to  early  failure. 

Problems  of  aelamination  growth  in  composites  have  recieved  attention 
in  recent  years. 

A  double  cantilever  beam  (D  C  B)  specimen  loaded  in  compression  was 
used  by  Chow  &  Ngan  (1)  to  study  the  effect  of  buckling  on  crack 
propagation.  Initial  geometric  imperfection  was  made  in  the  specimen  in 
order  to  obtain  crack  propagation  at  loads  below  the  buckling  load  of  the 
perfect  column.  Chai  (2)  carried  out  an  experimental  program  to  determine 

the  damage  propagation  mechanisms  of  a  composite  panel  (graphite/epoxy) 

•  • 

under  simultaneous  compressive  in-plane  loading  and  low-velocity  transverse 
impact.  Also,  he  presented  an  analytical  (one-dimensional)  model  for 
delamination  buckling  and  damage  growth.  For  a  homogeneous  isotropic  plate 
under  a  compressive  displacement-controlled  loading,  he  obtained  a  closed 
form  expression  for  the  energy  release  rate  for  the  case  of  a  thin  film 
model  (the  buckling  of  the  delaminated  layer  assumed  to  produce  a 
negligible  bending  deformation  in  the  main  body  of  the  laminate).  A  more 
general  model,  where  the  delamination  thickness  was  of  the  same  order  of 
magnitude  as  the  other  thicknesses  involved,  was  also  considered. 

Whitcomb  (3)  presented  an  approximate  analysis  for  a  symmetric 
through-width  delamination  of  a  unidirectional  graphite/epoxy  bonded  to  an 
aluminum  bar.  The  buckled  region  was  assumed  to  have  zero  slope  at  both 
ends.  For  the  case  of  short  and  thick  delaminated  regions,  the  rotations 
at  the  ends  cannot  be  neglected,  hence  he  accounted  for  the  rotations  by 
considering  a  modified  delaminatin  length  (effective  length)  which  he 
obtained  through  a  finite  element  solution  of  the  postbuckling  problem. 
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A  finite  element  stress  analysis  was  developed  (4)  to  study  the 
postbuckling  behavior.  The  stress  distribution  and  strain  energy  release 
rate  were  calculated  for  various  delamination  lengths,  delamination  depths, 
applied  loads  and  lateral  deflections.  A  few  specimens,  consisting  of 
unidirectional  graphite/epoxy  bonded  to  aluminum,  were  fatigue-tested  to 
obtain  delamination  growth  data.  Yin  and  Wang  (5)  derived  a  simple 
expression  for  the  energy  release  rate  associated  with  the  growth  of  a 
general  one-dimensional  delamination.  First,  a  continuous  stress  field 
produced  by  a  certain  subsystem  of  forces  and  moments  was  eliminated,  then 
the  energy  release  rate  associated  with  the  remaining  system  of  forces  and 
moments  was  evaluated  by  means  of  the  path-independent  J-integral.  Also, 
they  (5)  studied  the  possibilit  of  cracks  branching  off  to  neighboring 
interlaminar  faces  and  the  condition  for  coalescing  of  multiple  parallel 
delaminations  into  a  single  layer. 

A  two-dimensional  model,  consisting  of  an  internal  rectangular 
delarrunation  with  initial  displacement  imperfection  in  the  thinner  portion, 
was  proposed  by  Konishi  (6).  A  test  program  was  carried  on  a  hybrid 
specimen,  boron/graphite/epoxy  covers,  and  graphite/epoxy  substructure. 
The  effect  of  static  and  cyclic  loading  on  delamination  growth  was  studied. 
Moreover,  the  effects  of  temperature  and  moisture  were  taken  into 
consideration. 

The  case  of  a  circular  homogeneous  isotropic  plate  with  a  concentric 
penny-shaped  delaminatin  subject  to  a  radial  compressive  loading  was 
addressed  in  refs.  (7~9).  Based  on  an  asymptotic  analysis  of  the 
postbuckling  behavior  of  the  delaminate  layers,  Bottega  and  Maewal  (8) 
considered  the  case  of  symmetric  split,  where  the  delaminatin  assumed  to 
exist  and  grow  in  the  mid-plane  of  the  laminate.  Considering 
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time-dependent  loading  and  taking  into  account  the  inertia  of  the 
delaminated  layers,  Bottega  and  Maewal  (9)  examined  the  dynamics  of  the 
growth  of  penny-shaped  delamination  in  a  circular  plate.  The  load  applied 
was  in  the  form  of  a  radial  compressive  load  and  a  distributed  transverse 
tensile  pulse.  Using  Hamilton's  principle  the  equations  of  motion  and  the 
corresponding  growth  laws  were  derived.  Two  models  were  studied  by  Yin  (7) 
concerning  the  delamination  buckling  of  a  penny-shaped  delamination,  the 
thin  film  model  and  the  symmetric  split  model.  The  energy  release  rate  was 
obtained  for  both  models  (7)  by  means  of  the  M-integral. 

Using  finite  element  and  Rayleigh-Ritz  methods  Shivakumar  and  Whitcomb 
(10)  predicted  the  buckling  behavior  of  an  elliptic  delamination  embedded 
near  the  surface  of  a  thick  quasi-isotropic  laminate.  The  laminate  was 
symmetric  made  of  graphite/epoxy  laminae.  The  results  obtained  by  using 
Rayleigh  Ritz  analysis  showed  good  agreement  with  those  obtained  by  using 
finite  element  methods,  except  for  highly  anisotropic  sublaminates. 

The  present  work  deals  primarily  with  the  question  of  delamination 
buckling  and  growth  and  how  the  presence  of  the  delamination  affects  the 
global  load  carrying  capacity  of  the  system.  A  one-dimensional  model  is 
presented  in  order  to  predict  (a)  delamination  buckling  loads  of  across  the 
width  delaminated,  axially  loaded  laminated  plate  and  (b)  the  ultimate 
load-carry ing  capability  of  above  geometry,  when  delamination  growth  takes 
place.  The  model  is  employed  to  predict  critical  loads  for  delaminated 
wide  columns  with  both  simply  supported  and  clamped  ends.  The  effects  of 
delamination  position,  size  and  thickness  on  the  critical  loads  are  studied 
in  detail  for  both  sets  of  boundary  conditions.  Results  are  obtained  for 
delaminated  isotropic  plates  and  for  plates  made  up  of  special  type  of 
symmetric  cross-ply  laminates.  The  results  reveal  that  for  certain 


geometries  the  buckling  load  can  serve  as  a  measure  of  the  load  carrying 
capacity  of  the  delaminated  configuration.  In  other  cases,  the  buckling 
load  is  very  small  and  delamination  growth  is  strong  possibility,  depending 
on  the  toughness  of  the  material.  Also,  the  present  model  can  be  used  to 
study  the  effect  of  the  presence  of  coupling  between  bending  and  stretching 
on  delamination  growth.  In  this  aspect,  a  delaminated  plate  in  the  form  of 
unsymmetric  cross-ply  laminate  is  studied. 

1.2  Effect  of  Holes  and  Foreign  Inclusions  on  Buckling 

The  effect  of  small  cutouts  of  various  shapes  on  the  response  of 
elastc  shells  has  received  great  attention.  Experimental  studies  have 
shown  that  holes  could  cause  a  severe  reduction  in  the  buckling  loads  of 
shells. 

Based  ona  two-dimensional  finite  difference  scheme,  Brogan  and  Almroth 

(11)  presented  a  nonlinear  analysis  for  cylindrical  shells  with 
rectangular  cutouts.  A  modified  Newton  method  wa3  used  (11)  to  reduce  the 
computational  time.  Also,  simple  experiments  were  performed  (11)  in  order 
to  verify  the  computational  scheme  and  satisfactory  agreement  between  the 
theory  and  test  were  obtained.  A  finite  element  analysis  of  a  cylindrical 
shell  with  reinforced  circular  opening  was  presented  by  Meller  and  Bushnell 

(12) .  The  cylindrical  shell  was  ring-stiffened  and  made  of  steel,  and 
axially  compressed.  The  effect  of  geometric  imperfection  was  included  (12). 
At  the  same  time,  a  similar  work  to  that  done  by  Mellar  was  carried  out  by 
Baker  and  Bennett  (13).  They  (13)  tested  a  ring-stiffened  cylindrical 
shell,  reinforced  with  circular  holes  under  the  action  of  unsymmetric  axial 
loading.  The  effect  of  hole  diameter  on  the  buckling  load  was  reported 

(13) -  Also,  they  (13)  took  into  consideration  when  the  penetrations 
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(holes)  were  cutting  no  ring  stiffeners,  and  cutting  one,  two  and  three 
ring-stiffeners. 

As  far  as  the  effect  of  small  holes  on  composite  structure  is 
concerned,  an  experimental  program  was  carried  out  by  Byers  (14)  on  a 
damaged  graphi te/polymer  laminate.  A  controlled  damage  in  the  form  of 
circular  holes  and  simulated  delamination  were  introduced  into  the 
composite  specimens  (14),  which  were  subjected  to  static  compression  and 
cyclic  compression  load.  Four  different  types  of  graphite/polymer  were 
tested  to  evaluate  the  damage  tolerance  of  each  type.  In  some  cases  the 
compression  strength  of  the  damaged  specimens  reduced  to  40%  of  the 
compression  strength  of  the  undamaged  specimens. 

Another  material  imperfection  is  the  rigid  inclusion  (small).  The 
effect  of  rigid  inclusions  on  the  stress  field  of  the  medium  in  the 
neighborhood  of  the  inclusion  has  received  (limited)  attention  In  the  past 
25  years  [15-19].  As  far  as  these  authors  know,  the  effect  of  rigid 
inclusions  on  the  buckling  characteristics  of  the  system  (that  contains  the 
rigid  inclusion)  has  not  been  studied. 

1.3  Future  Work 

The  problem  of  buckling  of  delaminated  cylindrical  shells  has  not 
received  the  deserved  attention.  Very  few  investigations  have  been  carried 
out  in  this  area.  Kulkarni  and  Frederic  (20)  used  a  "branched  integration" 
technique  to  solve  the  problem  of  buckling  of  two  layered  cylindrical 
shell,  partially  debonded,  subjected  to  axial  compression.  They  (20) 
considered  the  case  where  the  delamination  originates  at  the  boundary. 
Results  were  obtained  for  different  lengths  of  debonding  and  inner  to  outer 
layer  thickness  ratios.  A  significant  decrease  in  the  critical  load  was 
observed.  The  buckling  of  stiffened  circular  cylindrical  shells,  with  two 
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unbonded  orthotropic  layers,  was  reported  by  Jones  (21).  Jones  (21) 
assumed  that  the  layers  did  not  separate  during  buckling,  i.e.  the 
deformation  of  both  layers  were  assumed  to  be  the  same.  Also,  he  examined 
the  case  where  one  of  the  two  unbonded  orthotropic  layers  was 
circumferentially  cracked.  His  results  for  a  cylindrical  shell  made  of 
aluminum  with  ablative  outer  layer  subjected  to  hydrostatic  pressure  showed 
that  the  ablative  layer  had  to  be  increased  by  50?  in  thickness  in  the 
damaged  (unbonded)  cylindrical  shell  in  order  to  obtain  the  same  buckling 
loads  of  the  perfect  cylindrical  shell.  Troshin  (22)  examined  the  effect 
of  longitudinal  delamination  in  a  laminar  cylindrical  shell  in  the  critical 
external  pressure.  The  delamination  was  assumed  to  extend  over  the  entire 
length  of  the  shell.  He  studied  the  effect  of  the  length  and  position  of 
the  delamination  over  the  shell  thickness  on  the  critical  external 
pressure. 

It  is  intended  to  extend  the  one-dimensional  model  of  sectin  (I)  to  a 
two-dimensional  model  to  investigate  and  study  the  effect  of 
circumferential  delamination  on  the  buckling  loads  of  circular  cylindrical 
shells  subject  to  compressive  loading.  Other  delamination  shapes  and  loads 
will  also  be  considered. 


II.  MATHEMATICAL  FORMULATION 
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II. 1  Description  of  the  Problem 

A  one-dimensional  modelling  of  the  laminated  plate  is  employed.  This 
is  similar  to  the  one  used  in  Ref.  3. 

Delamination  exists  and  grows  (if  it  does)  along  a  plane  parallel  to 
the  reference  plane.  This  simplified  modelling  is  employed  in  order  to 
understand  better  and  pay  more  attention  to  the  physical  behavior,  rather 
than  get  involved  (and  lost)  in  a  more  complex  model  with  several 
additional  parameters  and  a  considerably  more  difficult  solution  scheme. 
In  accordance  with  this  philosophy,  it  is  also  assumed  that  the 
delamination  exists  before  the  compressive  loading  is  applied.  The 
location  and  length  (size)  of  the  delamination  i3  arbitrary  and  the  ends  of 
the  plate  are  either  hinged  or  clamped  (see  Fig.  1).  Note  that,  because  of 
the  one-dimensional  modelling,  the  sides  of  the  plate  are  free.  Finally, 
it  is  assumed  that  the  delamination  separates  the  column  (or  plate)  into 
four  regions,  and  each  one  of  these  regions  has  such  dimensions  (lengths 
and  thicknesses)  that  Euler-Bernoull i  beam  theory  is  applicable  (this 
assumption  may  also  be  removed  by  employing  more  accurate  theories  for 
small  length  to  thickness  ratios). 

The  three  axes  of  orthotropy  are  parallel  to  the  reference  frame  x,y,z. 
The  natural  plane  of  the  plate  lies  on  the  xy-plane.  The  plate  is  of  unit 
width  and  is  subjected  to  a  uniform  compressive  load,  P,  along  the 
supported  edges.  A  layer  of  uniform  but  arbitrary  thickness  h  and  of 
length  a  is  delaminated  (3ee  Fig.  1).  The  delamination  extends  across  the 
entire  width  of  the  plate  (one-dimensional  modelling).  Note  from  Fig.  1 
that  the  delamination  divides  the  plate  (column)  into  four  regions,  1-4. 
The  coordinate  system  is  such  that  x  is  measured  from  the  left  end. 
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Moreover,  14  and  w^  (i  -  1,2,3,^)  denote  the  in-plane  and  transverse 
displacements  of  material  poins  on  the  midplane  of  each  region  (each  part). 
For  example,  W3  is  the  transverse  displacement  of  material  points  on  the 
plane  h/2  units  from  the  top  surface  (region  3),  while  W2  is  the  transverse 
displacement  of  material  points  on  the  plane  H/2  units  from  the  bottom 
surface  (region  2).  Furthermore,  note  that  because  of  the  Euler-Bernoulli 
assumptions  w^(x)  characterizes  the  transverse  displacement  of  every 
material  point  in  region  i  and  position  x.  Similarly,  the  inplane 
displacement,  ui(x,zj),  of  every  material  point  is  given  by 

^•(X,Z()=  3  w,tx  (j) 

where  z^  is  measured  from  the  midsurface  of  each  region  and  the  comma 
denotes  differentiation  with  respect  to  the  index  that  follows. 

II. 2  Delamination  Buckling 

The  delaminated  plate  (wide  column)  is  assumed  to  be: 

(a)  homogeneous  and,  at  most,  orthotropic  or 

(b)  made  up  of  special  type  of  symmetric  cross-ply  laminates;  such 
that  the  prebuckling  state  is  a  membrane  one. 

The  necessary  equations,  including  kinematic  relations,  constitutive 
relations,  relations  between  loads  and  moments  on  one  side  and  kinematic 
parameters  on  the  other,  and  equilibrium  equations,  are  given  below: 
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U) 


V/  -  -(&*;  K,x/ +  /Vxy.  U^JX  (5) 
=  0  (6) 


(6) 

(7) 


where 


axial  forces 
Mi  bending  moments 
Vi  shearing  force 

Note  that:  for  our  problem  (unit  width  plate) 


^t*i  "  ff 


where 


i  e»=°  [pM<  S^P) 


c  =  '  *  y  m 

**!  l  (ft  -J&L  \  j  so  [Wijt  Co/umfl) 

n  yy  ' 

/1TTS  J  Qjj  W 

-Lj  k- i  J  i  1,2, 1,4 


V+2  3>.'-C,) 
1 


U°) 


4 


4 


where  ,  V-^  and  P^  can  be  expressed  in  terms  of  the  displacement 
gradients,  through  Eqs.  (2)  to  (7). 

II. 2.1  Primary  (Membrane)  State  Solution 

As  the  compressive  load  P  is  applied  (quasi)  statically,  the  plate 
remains  flat  and  a  primary  state  solution  is  characterized  by 

wf-o 

W  </=4x  +  £,-  ^ 


where  and  are  arbitrary  constants. 

Use  of  the  boundary  and  continuity  conditions  leads  to  complete 
knowledge  of  the  primary  state,  or 
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The  related  conditions  are 
In-Plane  Boundary  Conditions 


Transverse  Boundary  Conditions 


(a)  Simple  Supports 

at  X-0  : 
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(b)  Clamped  Supports 

at  X=<?  ; 

\y,  -o 

a»cf 

II 

X 

n 

> 

Kinematic  Continuity 

Conditions 
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Continuity  in  Moments  and  Forces 
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T  .2.3  Solution  of  Buckling  Equations 

The  solution  to  the  buckling  equations  can  be  written  as 
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Note  that,  since  the  admissible  variations  in  the  displacement 
components  can  be  made  as  small  as  one  wishes,  the  reference  surface  strain 
can  be  linearized,  or  (see  Eq.  (2)). 


(2g ) 


Thus,  the  expressions  for  the  internal  forces  and  moments  become  (see 


Eqs.  (3-5), 
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Moreover,  the  solution  to  the  buckling  equations,  Eqs.  (27),  requires 
knowledge  of  24  constants  ( C  ^ i  ,  Cj2>  ai j >  1  1 >2,3,4  and  j  =  1, 2,3*4). 
There  exist  24  boundary  and  continuity  conditions,  which  are  homogeneous  in 
u^,  w^  and  their  space-derivatives.  These  consist  of  two  in-plane  boundary 
conditions,  Eqs.  (23),  four  transverse  boundary  conditions,  either  Eqs. 
(24a)  or  Eqs.  (24b),  twelve  kinematic  continuity  conditions,  Eqs.  (25a)  and 
(25b),  and  six  continuity  conditions  in  moments  and  forces,  Eqs.  (26a)  and 
(26b).  Their  number  is  24  also. 

Use  of  the  boundary  conditions  yields  a  system  of  linear,  homogeneous, 

algebraic  equations  in  C^,  and  ai j  •  For  a  nontrivial  solution  to 

exist  the  determinant  of  the  coefficients  must  vanish.  The  determinant 

contains  geometric  parameters  and  the  applied  load  P,  because  of  the 

P 

presence  of  [see  Eqs.  (18)]. 


Before  proceeding  with  the  description  and  actual  solution  of  the 
determinant,  certain  non-dimensionalized  parameters  are  introduced.  These 
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Note:  for  the  considered  case,  where  all  four  regions  (plates)  are 

symmetric  w.r.t.  their  mid-reference  planes,  the  following  relations  holds: 
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Now,  we  can  write  the  following  expressions  for  the  kjj's 
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Moreover,  the  nondimensionalized  coefficients  a<  . 

XJ,  are  defined  as 
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In  terms  of  the  new  non-dimensionalized  quantities,  the  system  of 
homogeneous,  linear  algebraic  equations  which  govern  the  buckling  of  the 
delaminated  plate,  could  be  written  as: 


(17) 


where  [A]  is  a  twenty  by  twenty  matrix  and  { X }  a  one  by  twenty  (column) 
matrix. 

I^f  ~  l*"  fin  fin,  c \  au 

&ll  #11  (Zj  f  4^2  ^7  ^4  ^7t  ^7 1  ^ J  ]?* J 

For  a  nontrivial  solution  to  exist  for  the  system  of  Eqs.  (37),  the  M) 
determinant  of  the  coefficients  must  vanish  (see  Eq.  39).  The  lowest 
eigenvalue  (P)  gives  the  value  of  the  critical  load.  Note  that  the 
elements  of  the  determinant  in  Eq.  (39)  contain  the  geometric  parameters  h, 
l  and  a  as  well  as  the  material  parameters,  Exx,  Eyy,  Vxy,  the  effect  of 

stacking  sequence  included  in  Dxx; ,  load  parameter  P. 

/ 

1.3  Delamlnatlon  Growth 

Delamination  of  a  composite  laminate  reduces  the  over-all  stiffness 
and  thereby  lowers  the  buckling  load  of  the  laminate.  The  latter  may  or 
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may  not  be  a  useful  indication  of  the  load-carrying  capacity  of  the 
delaminated  plate.  In  some  cases  it  may  be  significantly  smaller  than  the 
elastic  collapse  load  in  the  postbuckling  regime,  especially  if  the 
laminate  contains  a  relatively  long  and  thin  delamination.  In  order  to 
determine  the  ultimate  load  capacity,  and  the  possibility  of  spreading  of 
the  damaged  (delaminated)  area  to  the  undamaged  area,  the  post-buckling 
solution  is  employed.  Delamination  growth  (if  it  occurs)  is  assumed  to 
occur  in  its  own  plane  and,  it  is  governed  by  the  Griffith  fracture 
criterion. 

If  growth  takes  place,  it  is  important  to  examine  whether  the  growth 
could  be  arrested  at  a  later  stage  or  not.  The  energy  release  rate  has 
been  used  to  determine  whether  the  delamination  growth  is  stable  or 
unstable.  The  path-independent  J-integral  (27)  has  been  used  to  obtain  the 
energy  release  rate. 

Also,  the  present  model  will  be  used  to  study  the  effect  of  the 
presence  of  coupling  between  bending  and  stretching  on  delamination  growth. 
In  this  aspect,  a  delaminated  plate  in  the  form  of  unsymmetric  crossply 
laminate  is  studied.  Note  that  for  such  geometry  there  is  coupling  between 
bending  and  stretching  regardless  of  the  level  of  the  applied  load,  and 
therefore  the  possibility  of  bifurcational  buckling  does  not  exist. 

III. 3.1  Mathematical  Formulation  of  Delamination  Growth 

In  deriving  the  equations  governing  the  delamination  growth,  the  case 
of  symmetric  delamination,  where  the  delamination  is  located  symmetrically 
w.r.t.  the  two  ends  of  the  plate,  will  be  considered.  Also,  a 
claraped-clampled  boundary  condition  is  assumed.  In  the  following  analysis, 
the  case  of  plate  strip  will  be  considered  (e y  -  0).  The  model  will  handle 
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both  cases;  the  case  for  which  the  primary  state  is  membrane  one  and  the 
case  where  there  is  a  coupling  between  bending  and  stretching. 

The  stress  resultants  and  bending  moments  are  given  by  (26); 


A  -8  Is 
B  -D  I* 


(4°) 


where  for  the  plate  strip 


u;  o  o) 

(XlT=(^  o  o} 

Note  that  for  both  the  orthotropic  plate  and  cross-ply  laminate 

An  =  =  ty}  =  £>  «  o 


The  equilibrium  equations  for  plate  (i)  reduces  to: 


A ^  =-  'P 
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where 
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Substitution  of  the  expression  for  NXx  and  MXx  into  the  equilibrii 


equations,  Eqs.  (43,44)  yields 
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To  obtain  the  in-plane  displacement  we  have 


For  each  part  of  the  four  parts  (plates)  of  the  delaminated  plate  the 


in-plane  displacement  and  transverse  displacement  are  given  by: 

l/-CxJ-2y  ty,x 
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For  symmetric  delamination  with  clamped-damped  boundary  conditions,  Eq. 


(24b)  could  be  simplified  to: 
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The  displacements  and  (i  -  1,2,3)  are  now  given  in  terms  of  only 
three  unknowns  A,  K2  and  K3,  since  Ki  is  given  in  terms  of  the  applied  load 
P.  To  determine  these  three  unknowns,  three  equations  are  required.  The 
required  equations  are  the  continuity  in  forces  and  moments  as  well  as 
continuity  in  in-plane  displacement. 


-  continuity  in  forces 


CSV 


-  continuity  in  moments  (see  Fig.  2) 


(fXJ 


Using  the  expressions  for  Mi ,  m2  4  M3  gives 
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where 


0  =  flX{  5  i/t  k,/t 

-  continuity  in  in-plane  displacement  (see  Fig.  2) 
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Using  Eqs .  (49  &  58)  resuls  in: 
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Using  Eq.  (55),  we  have 
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Solving  Eqs .  (53,  59,  &  60)  simultaneously,  one  can  obtain  the  values  of 
P1 ,  ?2  &  P3  in  terms  of  the  applied  load  P.  Also,  once  Pi,  ?2  and  p3  are 
known,  the  unknown  constant  A  could  be  obtained  and  hence,  the  displacement 
functions  are  known  for  every  level  of  the  applied  load  P.  Now,  the 
postbuckling  solution  is  known  and  the  next  step  is  to  establish  a 
delamination  growth  criterion. 

II. 3. 2  ENERGY  RELEASE  RATE 

In  the  present  work,  it  is  assumed  that  delamination  growth  (if  it 
occurs)  will  take  place  in  its  own  plane  and  is  governed  by  the  Griffith 
fracture  criterion.  It  is  also  important  to  find  out  whether  this  growth 
is  stable  or  unstable,  in  doing  so,  one  has  to  calculate  the  energy  release 
rate  (the  energy  required  to  produce  a  unit  of  new  delamination). 

In  order  to  calculate  the  energy  release  rate,  the  path-independent 
J-integral  is  used  (27).  The  use  of  the  J-integral  beside  being  helpful  in 
avoiding  the  calculation  of  stress  singularities  at  the  delamination  tip, 
is  useful  in  the  cases  where  plasticity  effects  are  not  negligible  (28). 

The  J-integral  is  defined  as: 


where  r  is  a  closed  contour  around  the  crack  tip. 
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but  £y  a  0 

G  =  J=ij*u*‘/?  (63J 

where  G  is  the  energy  release  rate 


where  ox^  is  the  stress  in  the  x-direction  in  lamina  number  k  belongs  to 
region  (i). 

The  stress  distribution  is  decompassed  into  two  subsystems,  as  in  Refs 
(5,29),  (a)  a  continuous  stress  field,  which  makes  no  constribution  to  the 
J-integral  and  (b)  a  subsystem  (2)  which  does.  The  subsystem  (2) 
contributions  are: 

Region  (3)  contributions 
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II. 3. 3  End  Shortening  and  Delamination  Opening 


As  we  will  see  later,  the  knowledge  of  the  end  shorting  and 
delamination  opening  will  help  in  some  cases  to  find  the  load  carrying 
capacity  of  the  system. 
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First  the  following  non-dimensional  parameters  are  introduced. 
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In  terms  of  the  above  nondimenslonal  parameters,  the  equations  governing 
the  behavior,  Eqs  (53,  59,  &  60)  and  growth  Eq.  (68)  of  the  delaminated 
plate  can  be  summarized  as: 


Kl  ix  Ki  *'*  K 


These  three  equations,  Eqs.  (78)— (80)  are  valid  for  any  plate  with 
across  the  width  (one-dimensional)  delamination,  provided  that  the 
assumption  of  cylindrical  bending  deformation  holds  (ey  ■  0). 
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III  RESULTS  AND  DISCUSSION 

The  Georgia  Tech  high  speed  digital  computer  CDC-Cyber  70,  Model  74-28 
was  employed  for  generating  results.  These  results  are  presented  both 
in  tabular  and  graphical  form. 

III.1  Buckling  Results 
III.  1.1  Orthotropic  Plate 

Results  for  buckling  loads  of  a  delaminated  orthotropic  plate  are 
generated. 

Table  1  shows  values  of  critical  loads,  p,  of  a  clamped  configuration 
with  a  symmetric  delamination  [see  Eq.  (30)],  for  several  values  of  the 
delamination  length  parameter,  a",  and  of  the  delamination  thickness 
parameter,  h. 

The  last  row  gives  the  sum  of  the  buckling  loads  of  the  two  parts,  if 
the  delaraination  extends  through  the  entire  length  a  »  1  .  This  sum,  p,  is 
a  measure  of  the  load  carrying  capacity  of  the  completely  damaged 
configuration 

Pc^/r,=  h +U-A)  (Si) 

These  results  are  also  shown  graphically*  in  Fig.  3.  It  is  clearly 
seen  from  these  results,  that  as  long  as  a  S  h,  buckling  load  is  not 
affected  appreciably  by  the  presence  of  the  delaminatior.  provided  that  the 
delamination  thickness  is  relatively  small  (h  <  .2).  For  the  same  condition 
of  a  £  h,  the  presence  of  the  delamination  becomes  more  pronounced  as  one 
approaches  a  raidthickness  delaraination  (h  -  .5).  Note  that  for  this 
extreme  case,  when  a  »  h  -  .5,  Pcr  ■  .69  and  Pcr.  increases  as  a  becomes 

smaller  and  smaller.  On  the  other  hand,  when  I  >  h  the  buckling  load  is 
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Table  1.  Buckling  Loads  for  a  Symmetric  Delamination-Clamped  Boundaries 


\  h 
ffNy 

.02 

.05 

.10 

.20 

.30 

.40 

.50 

.025 

.6399889 

.9999989 

.9999996 

.9999999 

.9999999 

.9999999 

•9999999 

.05 

.1599975 

.9949315 

.9999824 

.9999929 

.9999949 

.9999938 

.9999960 

.10 

.0399994 

.2499388 

.9799200 

.9997142 

.9998254 

.9998600 

.9996683 

.15 

.0177775 

.1110863 

.4434620 

•9965731 

.9985215 

.9988911 

.9989777 

.20 

.0099999 

.0624871 

.2495341 

.9264435 

.9923920 

•9950435 

.9955576 

.25 

.0063999 

.0399923 

.1597332 

.6255729 

.9661928 

.9635314 

.9859276 

.50 

.0044444 

.0277728 

.1109427 

.4371092 

.8582246 

.9543270 

.9638297 

.40 

.0025000 

.0152263 

.0624213 

.2470476 

.5313834 

.7883057 

.8480737 

.50 

.0016000 

.0099987 

.0399587 

.1585467 

.3469048 

.5675410 

.6895626 

.60 

.0011111 

.0069437 

.0277551 

.1103405 

.2455310 

.4123884 

.5411401 

.70 

.0008163 

.0051017 

.0203958 

.0812210 

.1803553 

.3111042 

.4309655 

.80 

.0006250 

.0039061 

.0156187 

.0622944 

.1390114 

.2428107 

.3514201 

•  90 

.0004938 

.0030865 

.0123432 

.0493028 

.1104955 

.1948719 

.2923016 

1.00 

.0004000 

.0025000 

.0100000 

.04000X 

.090X00 

.160X00 

.25000X 

rdam 

.9412000 

.6572000 

.7300x0 

. 5200000 

.3700X0 

.28XX0 

.250X00 

greatly  affected  and  it  drops  drastically,  especially  for  the  smaller 
values  of  the  delamination  thickness. 

A  comparison  of  the  buckling  loads  with  the  buckling  load  for  the 
completely  damaged  system  (last  number  in  each  column)  reveals  certain 
observations  which  suggest  related  conclusions.  When  a  £  h  the  values  for 
Pcr  are  higher  than  P^am*  This  suggest  that  for  these  geometries  the  load 
carrying  capacity  of  the  system  is  related  to  and  measured  by  the  buckling 
load  of  the  delaminated  configuration.  On  the  other  hand,  for  a  >  h 
buckling  occurs  at  low  values  of  the  applied  load.  Then  if  the  load  is 

further  increased  this  may  lead  to  delamination  growth  and  the  damaged  area 

•  • 

may  extend  along  the  enire  length  of  the  system.  This  of  course  depends  on 
the  fracture  toughness  of  the  material,  but  the  implication  is  that  Pdam  i3 
a  good  lower  bound  for  the  load  carrying  capacity  of  the  delaminated 
conf iguration.  The  problem  of  postbuckling  behavior,  delamination  growth 
and  accurate  estimation  of  the  load  carrying  capacity  of  a  delaminated 
configuration  is  the  subject  of  Part  III. 2  of  this  report. 

On  Table  2,  the  value  of  the  load  in  region  3  at  the  Instant  of 
buckling  is  given  for  the  same  range  of  parameters,  a  and  h,  as  on  Table  1. 
These  results  are  also  for  the  case  of  clamped  supports.  The  values  for 
this  region  3  load  are  nondimensionalized  with  respect  to  the  critical  load 
of  the  region  3  geometry,  as  if  its  ends  were  clamped.  This  is  done 
primarily  for  finding  the  range  of  parameters  a  and  h  for  which  thin  film 
analysis  holds  [2].  Clearly,  when  pj  is  close  to  unity,  thin  film  analysis 
is  applicable  and  thus  delamination  growth  can  be  treated  by  the  simpler 
analysis  [2]  developed  for  the  case  of  thin  film  behavior.  Thin  film 
behavior  implies  that  region  3  only  experiences  buckling  and  postbuckling 
deformations  while  the  remaining  of  the  plate  remains  undeforraed.  For 


Table  2.  Region  of  Applicability  of  Thin  Film  Analysis  for  a  Symmetric 
Delaminated,  Wide  Column  with  Clamped  Ends 


\h 

*  \ 

.02 

.05 

.10 

.20 

.30 

.40 

.50 

.025 

.9999827 

.2499997 

.0625000 

.0156250 

.0069444 

.0059062 

.0025000 

.05 

.9999045 

.9949313 

.2499956 

.0625000 

.0277776 

.0156249 

.0100000 

.10 

.9999855 

.9997555 

.9799200 

.2499286 

.1110917 

.0624912 

.0399947 

.15 

.9999863 

.9997769 

.9977894 

.560572$ 

.  496303 

.1404691 

.0899000 

.20 

.9999d71 

.9997932 

.9981364 

.9264439 

.4410631 

.2487609 

.1592892 

.25 

.9999801 

.9998079 

.9983327 

.9774573. 

.6709672 

.3841920 

.2464012 

.50 

.9999609 

.9998216 

.9984847 

.9834956 

.8502246 

.5360090 

•3469787 

.40 

.9999905 

.9998483 

.9987409 

.9081904 

.9446817 

.7803037 

.5427671 

.50 

.9999923 

.9998741 

.9989605 

.9909166 

.9636244 

.0867828 

.6095626 

.60  . 

.9999937 

.9998995 

.9991833 

.9930647 

.9741242 

.9270740 

.7792417 

.70 

.9999953 

.9999261 

.9993921 

.9949577 

.9819344 

.9527566 

.8446920 

.80 

.9999970 

.9999503 

.9995985 

.9967103 

.9885257 

.9712427 

.8996356 

•90 

.9999985 

.9999750 

.9997994 

.9983807 

•9944594 

.9865390 

.9502971 

1.00 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 
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small  values  of  h,  say  h  ^  0.02,  thin  film  is  applicable  for  all  values  of 
a.  As  h  Increases,  towards  its  maximum  value  of  0.5,  the  region  of 
applicability  of  thin  film  analysis  is  confined  towards  the  high  a-value 
region.  The  results  of  Table  2  are  presented  graphically  on  Fig.  4. 

Similary  results  are  presented  on  Tables  3  and  4  and  Figures  5  and  6 
respectively  for  the  case  of  simply  supported  ends.  On  Table  3.  P  is 
nondimensinalized  with  respect  to  the  perfect  geometry  critical  load  for 
simply  supported  ends.  The  values  corresponding  of  Pdam  (last  row),  are 
the  sum  of  the  contributions  of  both  parts  acting  independent  of  each  other 
and  both  being  simply  supported.  This  may  be  questionable,  because  a 
better  description  of  the  completely  damaged  plate  is  needed.  For  example, 
if,  when  the  delamination  extends  along  the  entire  length,  the  thin  part  is 
lost  (it  becomes  disconnected)  then  a  better  measure  for  pdaffl  would  be 
(1-h)3  rather  than  h3  +  (1-fi)3.  it  is  also  worth  mentinoing  here  that  the 
values  (P)  corresponding  to  a  ■  1  (next  to  the  last  row)  are  calculated' 
under  the  assumption  that  both  parts  (regions  3  and  2)  have  the  same  slope 
at  the  boundaries. 

The  trends,  observed  in  the  case  of  clamped  supports,  are  also 
observed  here  but  She  regions  seem  to  have  shifted.  For  instance,  on  Table 
1  it  was  observed  that,  as  long  as  a  S  h,  the  buckling  load  was  not 
affected  appreciably  by  the  presence  of  the  delamination,  for  relatively 
small  delamination  thicknesses.  Here  on  Table  3  this  statement  is  also 
true  but  a  £  h  must  be  replaced  by  I  S  2h,  and  so  on. 

So  far  results  have  been  presented  and  discussed  for  the  case  of 
symmetrically  located  delaminations  [ii  -  0.5(1 -a ) ] . 

Next,  results  are  presented  on  Tables  5,6,7  and  8  for  the  case  for 
which  the  delamination  is  asymmetrically  located. 
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Table  3.  Buckling  Loads  for  a  Symmetric  Delamination  Simply  Supported 
Boundaries 


B 

wm 

.05 

•  io 

.20 

.30 

.40 

.50 

.02  5 

•  Jj'jjyS  9 

X • JoOvUO 

1  •  Gv^OuCG 

1 . 000000 

1.000000 

1 .000000 

.05 

.'j599u7o 

.999997a 

.9999992 

.9999996 

.9999997 

.9999999 

•9999999 

•  ID 

.1599975 

•9929014 

.3999 to  5 

.9999051 

.9999897 

.9999914 

.9999919 

.15 

.0711099 

.4445054 

•9995550 

.9990786 

.9999198 

.9999542 

.9999379 

.20 

.0599994 

.2499286 

.9722709 

•9994209 

.999o506 

.9997168 

.9997368 

.25 

.0255996 

.1399562 

.6376529 

.9979552 

.9966652 

•9991293 

.9991891 

.50 

.0177775 

.1110815 

.4431942 

•3952354 

.9970534 

.9977360 

.9979614 

.40 

.0099996 

.0624636 

.2493640 

.9016443 

.9651533 

.9902270 

.9912249 

.50 

.0064000 

.0 399897 

.1596320 

.6177589 

.9402253 

.9686038 

.9729376 

.60 

.0044444 

.0277707 

.1108654 

.4530717 

.8149228 

.9198068 

.9343196 

.70 

.0032653 

.0204030 

.0614568 

.3195488 

.6484304 

.8350811 

.8703371 

.60 

.0025000 

.0156210 

.0623679 

.2449950 

.5117912 

.7263620 

.7867299 

.90 

.0019753 

.0123426 

.0492799 

.1938276 

.4105782 

.6176623 

.6966477 

1.00 

.0016000 

.0099975 

.0599176 

.1571446 

.5556235 

.522  8  852 

.6109859 

w 

•  941^jj0 

.8572000 

.7300000 

. 5200000 

.5700000 

.2800000 

.2500000 

Table  4  Region  of  Applicability  of  Thin  Film  Analysis  for  a  Symmetric 
Delaminated.  Wide  Column  with  Simply  Supported  Ends 


X 

.02 

.05 

.10 

.20 

.30 

.40 

.50 

•025 

.3906251 

.0625000 

.0156250 

.0039063 

.0017361 

.0009766 

.0006250 

-05 

.9999810 

.2499995 

.0625000 

.0156250 

.0069444 

.0039062 

.0025000 

.10. 

.9999834 

.9929014 

.2499913 

.0624991 

.0277775 

.0156249 

.0099999 

.15 

•9999856 

.9996871 

.5622497 

.1406079 

.0629500 

.0351539 

.0224986 

.20 

.9999838 

.9997154 

.9722769 

.2498572 

.1110723 

.0624824 

.0399895 

.25 

.9999839 

.9997260 

.9963014 

.3898177 

.1734176 

.0975712 

.0624495 

.30 

.9999840 

.9997337 

•9970970 

.5586949 

.2492658 

.1403147 

.0898165 

.40 

.9999841 

.9997383 

.9975360 

.9016460 

.4378460 

.2475568 

.1585960 

.50 

.9999842 

.9997416 

.9977000 

.9652485 

.6529328 

.3783609 

.2432344 

.60 

.9999843 

.9997439 

.9977883 

.9744113 

.6149228 

.5173914 

.3365550 

.70 

.9999643 

.9997459 

.9978456 

.9780056 

.8825858 

.6593590 

.4264652 

.00 

.9999843 

.9997466 

.9978860 

.9799799 

.9098510 

.7263620 

.5055072 

.90 

.9999843 

.9997483 

.9979178 

.9012521 

.9238010 

.7817289 

.5642846 

1.0Q 

.9999843 

.9997483 

.9979401 

.9821539 

.9322875 

.81700 50 

.6109859 

Fig.  6  Region  of  Applicability  of  Thin  Film  Analysis  for 
Symmetric  Delamination.  Wide  Column  with  Simply 
Supported  Ends . 


Table  5.  Effect  of  Delamination  Location  on  the  Buckling  Load 
(Clamped  Supports;  a  =  .1) 


.05 

HU 

mm 

.3 

•4 

o - 

.0 

.2499536 

•960b663 

.9991479 

.9967055 

.998080? 

.9977440 

.1 

.2499509 

.9670979 

.9956626 

.9922549 

.9679301 

.9656805 

•2 

.249946} 

.997 -;b2  4 

.9938936 

.9864996 

.9822217 

•9790165 

•  5 

.2499416 

•9666065 

.9958984 

.9923642 

.9862422 

.9661157 

.45 

.2499366 

.9799206' 

.9997142 

.999  62  54 

.9998600 

.99986S3 

Table  6.  Effect  of  Delamination  Location  on  the  Buckling  Load 
(Clamped  Supports;  a  *  .2) 


iffii 

.05 

.1 

.2 

.3 

.4 

.5 

n 

.0624932 

.2497486 

•9527570 

! 

I  -9557411 

.9247731 

.9051069 

.0624907 

.2496894 

I  .9479818 

.9142127 

.8652345 

.6405072 

n 

.0624608 

.2496161 

.9455077 

.9194764 

.8767561 

.8554b76 

•  3 

.0624876 

.2495566 

.9309424 

.9613436 

.9401658 

.9277658 

n 

.0624871 

.2495341 

.9264435 

.9925920 

.9950435 

•9955576 
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Table  7.  Effect  of  Delamination  Location  on  the  Buckling  Load 


(Simply  Supports;  a  -  .1) 


S3 

.0 

.1 

.2 

.3 

.45 

.05 

.998639 

.996576 

.994855 

.993624 

.992901 

.10 

.999407 

.999511 

.999679 

.999847 

.999965 

.20 

.998517 

.998790 

.999232 

.999674 

.999985 

.30 

.997760 

.997760 

.998583 

.999409 

.999990 

.40 

.995823 

.996598 

.997849 

.999106 

.999991 

.50 

.995113 

.996020 

.997483 

.998954 

.999992 

Table  8.  Effect  of  Delamination  Location  on  the  Buckling  Load 
(Simply  Supports;  a  -  .2) 


F4 

.0 

.1 

.2 

.3 

.4 

.05 

.249936 

.269933 

.269931 

.249929 

.249929 

.10 

.988205 

.981420 

.976383 

.973308 

.972277 

.20 

.989043 

.991885 

.995416 

.998312 

.999429 

.30 

.979666 

.985074 

.991836 

.997459 

.999651 

.40 

.969003 

.977209 

.987545 

.996274 

.999719 

.50 

.963667 

.973241 

.985342 

.995645 

.999737 
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Tables  5  and  6  present  results  for  the  case  of  clamped  supports,  in 
terms  of  buckling  loads  p  for  various  values  of  h,  and  l\  .  The  results  on 
Table  5  correspond  to  I  -  0.1  and  those  on  Table  6,  to  I  «  0.2.  It  is  seen 
from  Table  5  and  6  that,  as  long  as  h  £  a,  the  buckling  load  is  the 
smallest,  when  the  delamination  is  located  symmetrically,  with  respect  to 
the  midpoint  of  the  wide  column.  On  the  other  hand,  when  h  >  a  the 
buckling  load  is  the  smallest  when  the  delamination  is  locatd  in  such  a  way 
that  it  spans  the  quarter  point  of  the  wide  column. 

For  the  case  of  simply  supported  columns  the  results  of  Tables  7  and  8 
reveal  similar  trends.  For  h  <  a,  the  buckling  load  is  the  smallest  when 
the  delamination  is  located  symmetrically  with  respect  to  the  midpoint  of 
the  wide  column.  Finally,  for  h  i  a,  the  buckling  load  is  the  smallest, 
when  the  delamination  starts  from  the  end  of  the  wide  column. 

Note  that  in  some  cases,  for  h  >  a  ■  0.2  and  clamped  boundary 
conditions,  effect  of  delamination  on  the  buckling  load  is  considerable. 
From  Table  6  one  observes  that  for  h  -  0.5  pmin  "  0.841  (j^  -  0.1)  while 

Praax  ’  0-996  (T,  -  0.4). 

III. 1.3  Symmetric  Cross-ply  Plates 

Results  for  buckling  loads  of  delaminated  plates  made  up  of 
Graphs te/Epoxy  have  been  obtained. 

The  orthotropic  axes  are  alternately  oriented  at  angles  0°  and  90° 
with  the  structural  axes. 

The  elastic  constants  typical  of  this  material  are: 


where, 


£i/£t  =4-o  ,  G’Lr/eT-  S 


j 


*■25 


El  is  the  tensile  modulus  in  the  filament  direction  (30  x  10^  psi) 
Ef  is  the  modulus  in  the  transverse  direction  (.75  x  10^  psi) 
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Table  9.  Buckling  Loads  for  Symmetric  Cross-Ply  Laminates, 


[0°/90°/0°3^qtJ  with  Clamped  Ends 


1 

.3 

.4 

.5 

0.05 

•9999oa6 

•9999933 

.9999951 

.9999956 

.9999960 

0.10 

.9990517 

•9997470 

.9998327 

.9998625 

.9998705 

0.15 

.6294934 

.9971932 

.9985950 

.99b9176 

.9989960 

0.20 

.3544494 

.9605042 

•9928885 

•9951751 

.9956416 

0.25 

.2269353 

.6856117 

•9695678 

.9840387 

.9862131 

0.50 

.1576354 

.4804679 

.8756095 

.9560119 

.9646097 

0.40 

.0od7059 

.2718917 

. 5517876 

.7965493 

.8556681 

0.50 

.0567913 

.1745731 

.3610403 

.5774029 

.6943155 

0.60 

.0394509 

.1215207 

.2536700 

.4205562 

.5460455 

0.70 

.0289931 

.0894789 

.1879582 

.3176146 

.4354499 

0.80 

.0222044 

.0686417 

.1449236 

.246O324 

.3554099 

0.90 

.0175495 

.0545362 

.1152281 

.1991666 

.2968545 

1.00 

.0142190 

.0440911 

.0938781 

.1635797 

.2531960 
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Table  10.  Buckling  Loads  for  Symmetric  Cross-ply  Laminates, 


[ 0°/ 90^/0°] 1Q  ,  with  Clamped  Ends 


.1 

.2 

.3 

.4 

.5 

0.05 

.9999907 

.9999937 

.9999953 

.9999960 

.9999961 

0.10 

.9994218 

.999764S 

•9998374 

.9998647 

.9998721 

0.15 

.7516396 

.9974912 

.9986399 

.9989348 

.9990081 

0.20 

.4237520 

.9714079 

.9931818 

.9952590 

.9956966 

0.25 

.2713595 

.7242170 

.9714534 

.9843584 

.9863974 

0.30 

.1885120 

.5039376 

.8856250 

.9570613 

.9651120 

0.40 

.1060935 

.2882870 

.5650983 

.8017537 

.8606837 

0.50 

.0679288 

.1851633 

.3703492 

.5838284 

.6974178 

0.60 

.0471911 

.1289302 

.2603643 

.4259284 

.5492680 

0.70 

.0346837 

.0949420 

.1929842 

.3219117 

.4384057 

0.80 

.0265641 

.0728424 

.1486349 

.2515180 

.3580450 

0.90 

.0209962 

.0576684 

.1183607 

.2020089 

.2992045 

1.00 

.0170128 

.0468003 

.0964461 

.1659502 

.2553130 
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Table  11.  Buckling  Loads  for  Symmetric  Cross-ply  Laminates, 
[ 0o/90°/0°]^q^,  with  Clamped  Ends 


.1 

.2 

.5 

.4 

•  5 

.9999925 

•999994*-* 

.9999955 

.9999959 

.9999961 

0.10 

.9995946 

.9997816 

.9996421 

.9998666 

.9998756 

0.15 

.8850965 

.9977555 

.9936859 

.9989527 

.9990209 

0.20 

.5005655 

.9784464 

.9954780 

.9.955480 

.99575^6 

0.25 

.520524a 

.7654652 

.9752845 

.9046975 

.9865968 

0.50 

.2226972 

.5402540 

.8954554 

.9581666 

.9656552 

0.40 

.1255515 

.5064068 

.5796544 

.8072950 

.8631029 

0.50 

.0802672 

.1968816 

.5806192 

.5908525 

.7008155 

0.60 

.0557675 

.1571250 

.2677667 

.4510495 

.5528164 

0.70 

.0409399 

.1009955 

.1985476 

.5266621 

.4416697 

0.80 

.0515961 

.0774970 

.1551672 

.2555549 

.3609600 

0.90 

.0243171 

.0615616 

.1218522 

.2051585 

.3018072 

1.00 

.0201110 

.0490056 

.0992950 

.1685782 

.2576600 

G^t  is  the  shear  modulus  (.375  x  10g  psi) 
is  Poisson's  ratio 

Results  are  generated  for  laminates  with  the  stacking  sequences, 
[0°/90°/0°]iot»  [0°/9QO/0°] ^ qt .  The  effect  of  delamination  length  and 
thickness  on  those  buckling  loads  is  studied.  Note  that  [0°/90i°/0°]ioT 
means  that  there  is  a  stacking  which  consists  of  one  thickness  of  0°,  i 
thicknesses  of  90°  and  one  thickness  of  0°,  and  repeated  ten  times. 

Results  obtained  for  symmetric  cross-ply  delaminated  plates  are 
similar  to  those  for  the  orthotropic  plate. 

Table  9  shows  the  values  of  the  buckling  loads,  F,  of  a  clamped 
symmetric  cross-ply  with  stacking  sequence  [0°/900/10°]iot»  for  several 
values  of  the  delamination  length  parameter  I,  and  of  the  delamination 
thickness  prameter,  h.  The  same  results  are  shown  graphically  in  Fig.  7. 
The  results  show  that  for  a  relatively  thin  delamlnatin  (h  S  .2),  the 
presence  of  delamination  has  a  negligible  effect  on  the  buckling  load  of 
the  delaminated  plate,  as  long  as  a  £  h.  On  the  other  hand,  for  the  case 
where  a  >  h  the  value  of  the  buckling  loads  is  greatly  affected  by  the 
presence  of  delamination,  especially  for  plates  with  thin  delamination. 

Tables  10  and  11  show  similar  results  for  the  same  configuration  but 
with  stacking  sequence  [0°/9Q^/0°]  and  [0°/90^/0°]  respectively.  The 
results  of  Tables  9,  10  and  11  are  compared  graphically  in  Fig.  8,  with  the 
results  obtained  previously  for  the  orthotropic  plate  (Table  1)  for 
delamlnatin  thickness  h  -.1. 

It  is  clear  from  Fig  8.  that  as  the  thickness  of  the  90°  layers 
increases  w.r.t.  the  thickness  of  the  0°  layers  the  value  of  the  buckling 
load  parameter  becomes  larger  for  the  same  delamination  length  and 
thickness.  On  the  other  hand,  we  have  to  realize  that  the  value  of  the 
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Table  12.  Buckling  Loads  for  Symmetric  Cross-ply  Laminates, 
[  with  Simply  Supports  Ends. 
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Fig.  10  Effect  of  Symmetric  Delamination  on  Buckling  Loads  of 
Cross-Ply  Laminates  with  Simply  Supported  Ends. 
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buckling  load  parameter  for  each  ca3e  is  normalized  w.r.t.  the  perfect 
configuration. 


The  effect  of  relative  value  of  transverse  modulus  (E22)  to  the 
tensile  modulus  (E^i)  is  studied.  Fig.  9  illustrates  this  effect  for  the 
case  of  a  symmetric,  delaminated,  cross-ply  plate  with  clamped-clamped 
boundary  condition  for  a  delamination  thickness  parameter  h  *  .1.  The 
considered  cross-ply  plate  has  stacking  sequence  [0/90/0]iqt*  The  results 
show  that  when  the  0°  layers  and  90°  layers  are  interchanged,  the  buckling 
load  parameter  drops  drastically  even  for  a  <  h.  A  similar  result  has  been 
obtained  for  a  [0°/90°/0°]  ^  g-p  cross-ply  laminate  with  simply  supported 
boundary  conditions.  These  results  are  presented  on  Table  12  and  in  Fig. 
10.  The  results  of  this  case  are  similar  to  those  of  the  orthotropic  plate 
in  general,  but  it  is  noticed  that  for  relatively  thin  delamination  (h  - 
.1)  the  buckling  load  parameter  for  the  cross-ply  laminate  is  much  higher 
than  that  of  the  orthotropic  plate.  This  difference  becomes  smaller  as  the 
delamination  thickness  increases. 

III.  1.3  Conclusions 

A  simple  model  has  been  employed  to  study  delamination  buckling  and 
the  effect  of  location,  size,  and  thickness  of  delamination  on  the  buckling 
load.  The  simplicity  of  the  model  limits  the  applicability  of  the  results 
to  the  case  for  which  each  region  of  the  four  parts  of  the  plate  is 
symmetric  w.r.t  its  reference  plane. 

For  these  geometries  (laminates)  the  results  serve  a  different  purpose. 
They  can  be  employed  to  relate  the  ultimate  load  carrying  capacity  of  the 
completely  (throughout)  damaged  laminate  to  the  allowable  for  the 


delaminated  laminate. 


A  good  application  of  the  present  analysis  and  results  is  in  the 
analysis  of  adhesive  bonding  of  two  similar  materials  [1], 

III. 2  Delamination  Growth  Results 

First  of  all,  before  presenting  any  results,  we  are  goint  to 
illustrate  the  growth  mechanism  of  the  delaminated  laminate,  let  us 
consider  Fig.  11.  Consider  a  material  with  critical  fracture  toughness  75* 
and  initial  delamination  length  a0.  The  load  is  increased  quasistatically 
to  pQ,  where  the  critical  value  of  the  energy  release  rate,  G*,  is  reached. 
As  the  delamination  length  extends  by  an  amount  Aa,  while  the  load  is  kept 

constant,  the  figure  shows  that  the  available  energy  release  rate 

•  • 

corresponding  to  the  new  delamination  length  a  ♦  Aa,  exceeds  the  critical 
value  G*,  hence  growth  will  take  place.  On  the  other  hand,  if  the  initial 
delamination  length  is  ai  ,  in  this  case  an  increase  in  the  delamination 
length  A31 ,  will  result  in  a  decrease  in  the  energy  release  rate, 
consequently  there  will  be  no  delamination  growth  as  long  as  the  load  does 
not  changed  (increased). 

III. 2.1  Orthotropic  Plate: 

Figures  12  to  14  contain  three  sets  of  curves  showing  the  relations 
between  the  non-diraensionalized  eneryg-release  rate,  G  -  L1*G/(Et5),  and  the 
normalized  delamination  length  a  »  a/L,  under  fixed  axial  load  and 
delamination  thickness.  The  three  figures  correspond  respectively  to  the 
cases  h  »  0.2,  0.1  and  0.02.  These  values  of  the  normalized  thickness  are 
thought  to  be  representative  of  relatively  thick,  a  relatively  thin,  and  an 
extremely  thin  delamination.  Each  curve  in  these  figures  refers  to  a  fixed 
value  of  the  normalized  axial  load  P,  and  each  curve  is  obtained  from  a 
one-parameter  family  of  numerical  solutions  of  Eq.  (75).  In  the  case  of 
relatively  long  delaminations  (a/h  >  1),  G  varies  over  a  wide  range  of 
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values  as  a  changes.  Hence  the  curves  are  plotted  in  a  semi-logrithmic 
diagram. 

In  a  quasi-static  process  of  delamination  growth,  the  energy-release 
rate  G  maintains  a  constant  critical  value  G*  (the  non-dimensionalized 
fracture  toughness). 

Figures  12-14  show  that  if  the  fractue  toughness  is  relatively  small, 
say  G*  <  10“2  (for  very  thin  delamination  h  -  .02)  and  G*  <  10-1  (for  h  £ 
.1),  then  the  delamination  growth  is  unstable,  in  general,  although  it  may 
become  stable  shortly  before  the  state  of  complete  delamination  (a  -  1  )  . 
Hence,  if  G*  is  small,  delamination  growth  under  a  constant  axial  load?  (a 
dead  load)  is  generally  a  catastrophic  process.  When  a  delaminated  plate 
buckles  under  an  increasing  axial  load  P,  the  postbuckling  solution  follows 
a  vertical  path  on  the  G-a  curves  (a  -  constant)  until  the  curve  intersects 
the  horizontal  path  G  -  <5*.  The  value  of  P  at  the  intersecting  point  is 
the  ultimate  axial  load  capacity  of  the  plate.  Afterwards  delamination 
growth  starts  and  proceeds  catastrophically  under  a  constant  axial  load. 

If  the  length  of  existing  delamination  is  relatively  short  (I  <  h), 
then  G  never  attains  the  critical  value  G*  and,  consequently,  delamination 
growth  does  not  occur.  For  such  plates,  the  ultimate  axial  load  capacity 
is  not  governed  by  delamination  growth,  but  is  determined  by  its  elastic 
postbuckling  behavior.  The  critical  buckling  load  is  a  lower  bound  of  , 
and  a  close  estimate  for,  the  ultimate  axial  load  capacity. 

If  the  fracture  toughness  is  relatively  large,  say  G#>  .3  (for  very 
thin  delaminatin,  h  *  .02)  and  G*  >  1 .  (for  h  £  .1),  the  figures  show  that 
all  curves  lie  above  the  curve  corresponding  to  Pdam  ^dam  “  h3  +  (1-h)3j). 
Consequently,  if  h  is  small  and  G  is  large  (as  given  above),  the  critical 
(buckling)  load  of  a  completely  delaminated  plate  is  a  close  lower  bound  of 
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the  ultimate  load  carrying  capacity  of  the  laminate.  If  the  fracture 
toughness  lies  between  these  two  extreme  cases,  rel  tively  high  and 
relatively  low,  then  the  delamination  growth  and  the  load  carrying  capacity 
will  depend  upon  the  value  of  the  fracture  toughness  ”5*  as  well  as  on  the 
initial  delamination  length  parameter.  The  above  statement  is  true  as  long 
as  the  deflection  and  the  end  shortening  are  both  bounded  (see  Figures  18 
and  19).  It  is  clear  from  the  obtained  results  that  for  the  3ame  level  of 
the  applied  load  and  for  the  same  delamination  length,  the  energy  release 
rate  increases  as  the  delamination  thickness  increases,  i.e.  for  the  same 
material  the  possibility  of  delamination  growth  increases  as  the  thickness 
of  the  delamination  increases. 

III. 2. 2  Symmetric  Cross-Ply  Plate 

The  energy  release  rate  is  obtained  for  a  delaminated  plate  made  up  of 
Graphite/Epoxy  with  same  material  constants  as  given  in  the  results  of  the 
buckling  problem.  Among  the  different  stacking  sequence  configurations, 
which  we  deal  with  in  the  buckling  problem,  we  are  going  to  consider  the 
case  of  the  stacking  sequence  [0°/90o/0°]iqT*  Fig.  15,  shows  the  results 
for  the  cross-ply  laminate  with  delaraination  thickness  h  -  .2.  The  general 
shape  of  the  obtained  curves  is  similar  to  that  of  the  orthotropic  plate. 
The  only  difference  between  the  cross-ply  results  and  the  orthotropic 
results  is  that  for  the  same  delamination  length  and  delamination  thickness 
and  some  applied  load,  the  energy  release  rate  for  the  cross-ply  laminate 
13  smaller  (not  much)  than  that  of  the  orthotropic  plate.  Note  that,  when 
we  say  same  level  of  the  applied  load,  we  have  to  keep  in  mind  that  this 
load  is  normalized  w.r.t.  the  perfect  geometry  plate  in  each  case. 
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19  Deflection  of  mid-point  of  delaminated  orthotropic  plate. 
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III. 2. 3  Unsymmetric  Cross-ply 

Since  for  laminates,  in  general,  there  exists  a  coupling  between 
bending  ane  extension,  the  present  model  is  used  to  study  the  effect  of  the 
presence  of  coupling  between  bending  and  stretching  on  delamination  growth. 
In  this  aspect,  a  delaminated  plate  in  the  form  of  unsymmetric  cross-ply 
laminate  is  studied.  Vote  that  for  such  geometry,  there  is  coupling 
between  bending  and  stretching  regardless  of  the  level  of  the  applied  load, 
and  therefore  the  possibility  of  bif urcational  buckling  does  not  exist. 
The  unsymmetric  cross-ply,  [90°/0°]iot  is  made  up  of  Graphite/Epoxy.  Figs. 
16  &  17  show  the  energy  release  rate  vs  the  delamination  length  for  various 
levels  of  the  applied  loads  for  delamination  thickness  h  -  .1  and  h  -  .2, 
respectively.  Unlike  the  orthotropic  case,  Fig.  16  shows  that  even  for  a  < 
h  there  exist  an  energy  release  rate,  and  the  curves  suggest  unstable 
growth  in  this  range.  The  results  show  that  for  the  same  parameters,  the 
energy  release  rate  for  the  unsymmetric  cross-ply  is  less  than  that  of  the 
orthotropic  plate.  Fig.  17  gives  the  energy  release  rate  for  the 
unsymmetric  cross-ply  with  delaraination  thickness  h  «  .2,  the  results  are 
similar  to  that  obtained  for  the  orthotropic  and  (symmetric)  cross-ply 
plate. 

REFERENCES 

1.  Chow,  C.  L.  and  Vgan,  K.  M.  "Method  of  Fracture  Toughness  Evaluation  of 
Adhesive  Joints"  J.  Strain  Analysis,  Vol.  15,  No.  2,  1980  pp.  97-101. 

2.  Chai,  H.  Ph.D.  Thesis,  "The  Growth  of  Impact  Damage  in  compressively 
loaded  laminates",  California  Institute  of  Technology,  Pasadena, 
California,  1982. 

3.  Whitcomb,  J.  D.  "Approximate  Analysis  of  Postbucling  Through-Width 
Delamination",  NASA  Technical  Memorandum  TM-83147  (1981). 

4.  Whitcomb,  J.  D.  "Finite  Element  Analysis  of  Instability  Related  Delami¬ 
nation  Growth",  J.  Comp.  Materials,  Vol.  15,  1981,  pp.  403-426. 

5.  Tin,  W.  L. ,  and  Wang,  J.  T.  S.  "Postbuckling  Growth  of  a  One-Dimen¬ 
sional  Delamination.  Part- I  Evaluation  of  the  Energy  Release  Rate", 

J.  Appl.  Mechanics,  Submitted  for  Publication,  1983. 

72 


% 


« 


6.  Konishi ,  D.  Y.  and  Johnston,  W.  R.  "Fatigue  Effects  of  Delamination  and 
Strength  Degradation  in  Graphite/Epoxy  Laminates",  Composite  Materials: 
Testing  and  Design  (fifth  Conference),  ASTM  STP  674,  S.  W.  Tsai,  Ed. 
1979,  pp.  597-619. 

7.  Yin,  W.L.  "Axisymmetric  Buckling  and  Growth  of  a  Circulr  Delamination 
in  a  Compressed  Laminate",  Int.  J.  Solids  and  Structures,  Submitted  for 
Publication,  1983. 

8.  Bottega,  W.  J.  and  Maewal ,  A.  "Delamination  Buckling  and  Growth  in 
Laminates",  J.  Appl.  Mechanics,  Vol.  50,  No.  1,  1983,  pp.  184-189. 

9.  Bottega,  W.  J.  and  Maewal,  A.  "Dynamics  of  Delamination  Buckling", 
Proceedings  of  AIAA/ASME/ASCE/AHS  24th  SDM  Conference,  Paper  83-0873 
Lake  Tahoe,  Nevada  May  1983. 

10.  Shivakumar,  K.  N.  and  Whitcomb,  J.  D.  "Buckling  of  a  Sublaminate  in  a 
Quasi-Isotropic  Composite  Laminate"  NASA,  TM-85755  (1984). 

11.  Brogan,  F.  A.  and  Almroth,  B.  0.  "Buckling  of  Cylinders  with  Cutout3", 
AIAA  J.,  Vol.  8,  1970,  pp.  236-260. 

12.  Meller,  E.  and  Bshnell,  D.  "Buckling  of  Steel  Containment  Shells" 
Lockheed  Palo  Alto  Research  Labs,  CA,  CR-7863  (Dec.  1982). 

13.  Baker,  w.  and  Bennet,  J.  "Buckling  Investigation  of  Ring-Stiffened 
Cylindrical  Shells  with  Reinforced  Openings  under  Unsymmetrical  Axial 
Loads"  Los  Alamos  National  Laboratory  report  LA-9646-MS  (Feb.  1983). 

14.  Byers,  B.  A.  "Behavior  of  Damaged  Graphite/Epoxy  Laminates  under 
Compression  Loading",  NASA  CR-159293  (August  1980). 

15.  Eshelby,  J.  D.  "The  Determination  of  the  Elastic  Field  of  an  Ellip¬ 
soidal  Inclusion,  and  Related  Problems"  Proceedings  Roy.  Soc.,  Series 
A,  Vol  241,  London,  1957,  pp.  376-396. 

16.  Walpole,  L.  J.,  "The  Elastic  Field  of  an  Inclusion  in  an  Anisotropic 
Medium",  Proceedings  Roy.  Co.,  Series  A,  Vol.  300,  London  1967,  pp. 
270-289. 

17.  Van  Dyke,  P.,  "Stress  in  a  Cylindrical  Shell  with  A  Rigid  Inclusion" 

AIAA  J.,  Vol.  5,  No.  1,  Jan.  1967,  pp.  125-137. 

18.  Kinoshita,  N.,  and  Mura,  T.,  "Elastic  Fields  of  Inclusions  in  Aniso¬ 
tropic  Media"  Physlca  Status  Solldl,  Vol.  (a)5,  1971,  pp.  759-768. 

19.  Takao,  Y.,  Taya,  M.,  and  Chou,  T.  W.,  "Stress  Field  Due  to  a  Cylin¬ 
drical  Inclusion  with  Constant  Axial  Elgenstrain  in  an  Infinite  Elastic 
Body"  J.  Appl.  Mech.,  Vol.  48,  No.  4,  1981,  pp.  853-858. 

20.  Kulkarni,  S.  V.  and  Frederick,  D.  "Buckling  of  Partially  Debonded 
Layered  Cylindrical  Shells",  Proceedings  of  AIAA/ASME/SAE  14th  SDM 
Conference,  Williamsburg,  Virginia,  March  1973. 


73 


21.  Jones,  R.  M.  "Buckling  of  Stiffened  Two-Layered  Shells  of  Revolution 
with  a  Circumferentially  Cracked  Unbonded  Layer",  AIAA  J.  Vol.  7,  1969, 
pp.  1511-1517. 

22.  Troshin,  V.  P.,  "Effect  of  Longitudinal  Delamination  in  a  Laminar 
Cylindrical  Shell  on  the  Critical  External  Pressure",  J.  of  Comp. 
Materials,  Vol.  17,  No.  5,  1983,  pp.  563“567. 

23.  Bellman,  R.,  Perturbation  Techniques  in  Mathematics,  Physics  and 
Engineering,  Holt,  Rinehart  and  Winston,  New  York,  New  York,  1969. 

24.  Sewell,  M.  J.,  "The  Static  Perturbation  Technique  in  Buckling  Problems" 
J.  Mech.  and  Phys.  of  Solids,  Vol.  13,  1965,  pp.  247-254. 

25.  Simitses,  G.  J.,  Elastic  Stability  of  Structures,  Prentice-Hall  Inc., 
Englewood  Cliffs,  N.J.,  1976  (Ch.  7). 

26.  Calcote,  L.  R.  "The  Analysis  of  Laminated  Composite  Structures", 

Van  Nostrand  Reinhold  Company,  New  York,  1969  (Ch.  6). 

27.  Rice,  J.  R.  "A  Path  Independent  Integral  and  the  Approximate  Analysis 
of  Strain  Concentration  by  Notches  and  Cracks",  J.  Appl.  Mech.  1968 
pp.  379-386. 

28.  Brock,  D.,  Elementary  Engineering  Fracture  Mechanics,  Martinue  Nijhoff 
Publishers,  1982  (Ch.  5). 

29.  Whitcomb,  T.  D.  "Strain  Energy  Release  Rate  Analysis  of  Cyclic  Delami 
nation  Growth  in  Compressively  Loaded  Laminates"  NASA  Technical 
memorandum  84598  (Jan.,  1983). 


74 


